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Abstract 

We study the spectra of the bottomonium and B c states within perturbative QCD up 
to order a^. The O(Aqcd) renormalon cancellation between the static potential and 
the pole mass is performed in the e-expansion scheme. We extend our previous analysis 
by including the (dominant) effects of non-zero charm-quark mass in loops up to the 
next-to-leading non- vanishing order s 3 . We fix the 6-quark MS mass = m^ IS (m^ IS ) 
on T(15) and compute the higher levels. The effect of the charm mass decreases mj, 
by about 11 MeV and increases the n = 2 and n = 3 levels by about 70-100 MeV 
and 240-280 MeV, respectively. We provide an extensive quantitative analysis. The 
size of non-perturbative and higher order contributions is discussed by comparing the 
obtained predictions with the experimental data. An agreement of the perturbative 
predictions and the experimental data depends crucially on the precise value (inside 
the present error) of a s (M z ). We obtain mf s (m^ s ) = 4190 ± 20 ± 25 ± 3 MeV. 



1 Introduction 



Traditional methods for theoretical investigations of heavy-quarkonium spectra have been 
based on various phenomenological potential models in the frame of non-relativistic quan- 
tum mechanics. These phenomenological potentials were tuned to reproduce the observed 
quarkonium spectra (and some other physical observables) and have been (in most cases) 
successful in predicting physical observables of the quarkonia, such as leptonic widths and 
transition rates among different levels, besides reproducing the energy levels. They have 
elucidated various properties of the bound states and essentially established, from consis- 
tency, the non-relativistic nature of the quarks inside quarkonia. For some reviews we refer 
to [|IJ. An apparent deficit of these approaches is, however, the difficulty in relating the 
phenomenological parameters to the fundamental parameters of QCD. 

Recently we have reported in new results on the spectroscopy of heavy quarkonia 
(charmonium, bottomonium and B c ) computed in perturbative QCD up to order af, once 
the C(Aqcd) renormalon cancellation between the static potential and the pole mass [[| Q 
has been implemented. In order to realize this we have adopted the so-called e-expansion 
scheme ||. (In the present work we will use the same scheme.) The major results of that 
analysis have been as follows. (1) Once the cancellation of the leading renormalon has been 
incorporated, the perturbative series turns out to be convergent and to reproduce reasonably 
well the gross structure of the bottomonium spectrum at least up to some of the n = 3 levels. 
(2) The constraints on non-perturbative and higher-order contributions to the bottomonium 
spectrum, set by the comparison of the calculation with the experimental data, indicate that 
these are smaller than usually believed. 

In this paper we improve the analysis of || by including the effects of the non-zero 
charm-quark mass loops on the spectra of bottomonium and B c . The typical scales of the 
bottomonium states as well as of the B c (as obtained in 0) are either close to or less than 
the charm-quark mass. Hence, charm-mass effects are expected and actually turn out to be 
numerically not negligible. These effects have been computed for the pole-mass-MS-mass 
relation at leading (non-vanishing) order in and at the next-to-leading order in the limit 
of small charm mass in |7]]. In the binding-energy expansion they have been calculated at 
leading (non- vanishing) order in || and at the next-to-leading order, only for the T(1S) 
mass, in 0. In the present work we implement and discuss the above results and derive 
some new formulas valid for the inclusion of charm- mass effects at the next-to- leading order 
for all quarkonium excited states. This will allow us to fully take into account charm-mass 
effects in the spectra of bottomonium and B c up to order e 3 in the e-expansion. Finally, we 
re-examine whether, by including these effects, the above conclusions (1) and (2) of still 
hold. 

The paper is organized as follows. In Sec. |2| we review the theoretical framework and 
the results of Ref. M. In Sec. |3| we present all the formulas needed for the inclusion of the 
non-zero charm-mass effects up to order e 3 in the bottomonium and B c energy levels. In 
Sec. |] we perform the numerical analyses on the bottomonium and B c spectra that include 
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the above effects. In Sec. || we give some conclusions. 



2 Quarkonium Spectrum with Massless Quark Loops 

The heavy-quarkonium dynamics is characterized by at least two dynamically generated 
scales: the (soft) scale mv, of the order of the momentum transfer and the (ultrasoft) scale 
mv 2 of the order of the energy of the quark in the bound state (v being the typical size of 
the heavy-quark velocity in a quarkonium state). The magnitude and nature of the non- 
perturbative corrections depend on the relative size of Aqcd with respect to the dynamically 
generated scales. Different scenarios are then possible. These have been systematically 
studied in [[|, where the effective field theories corresponding to them, generically denoted 
as potential non-relativistic QCD [|Kj, have been constructed. In |2| we have discussed our 
work in relation to these scenarios. That discussion applies also here. In order to facilitate 
the reading, let us summarize the main points. We assume to be in a situation where 
mv 3> Aqcd and the system to be mainly Coulombic. In this scenario non-perturbative 
corrections may be carried by local or non-local gluon condensates. Due to the fact that, in 
both direct numerical evaluation turns out to be quite uncertain, we chose not to 

correct explicitly our perturbative formulas with them (as done, for instance, in [|ll|, |12|), 
but rather to infer their size, as well as that of the neglected higher-order perturbative 
contributions from the eventual disagreement (and uncertainties) of the perturbative results 
from the experimental data. Among the neglected higher-order perturbative contributions we 
also include the In a s contributions, which are known from |13| , but whose dominance over 
other O(o.g) corrections is unclear. Therefore, using the effective field theory language of ||, 
we will consider in the following heavy quarkonium as described by potential non-relativistic 
QCD in the perturbative regime up to order a*. This will correspond to just considering the 
quarkonium spectrum up to that order as obtained from standard perturbative QCD. 

2.1 Perturbative expansions 

In this section, in order to set up some basic formulas, we review the energy (mass) of a 
quarkonium state X made by two heavy quarks in the theory with n\ massless quarks only. 
Namely, we do not take into account, for the moment, the effects of non-zero masses in the 
light-quark loops, neither in the binding energy nor in the pole-mass expansion. These will 
be considered in Sec. [3|. 

We express the quarkonium mass as a series expansion in the MS coupling constant 
a s ni ^ (^) defined in the theory with ni massless quarks only. Consider the quantum-mechanical 
Hamiltonian 



H = mi )PO i e + m 2 , P oic + H + H 1 + H 2 + ■ ■ 
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whereQ 



p 2 4a^)(/i) 

H ° = 2^ — ' (2) 

Hi = (^)U^i + a^V (3) 



3 r V 47r 



(^M) 2 (u 2 + £) + 2(/^ + ^ofty + 

r3/ . (1 -X)^)^) / 2 1 



(l + 3a;)p 4 7t(1 + o;)q;W)(//) r3/ ^ (1 - x) aj'")^) 



32 m 3 , 3 m 2 . 

(l-x)o^>(/i) fg 2 _ 3 (S^ _ 47r _ 3)(5 3 (f) l _ ^ W 'W 

6 m 2 I r 3 r 5 3 J m r r 2 

(3-x)«^(/x) T g , a^)(//) / 1 1 



6m 2 H 3H \ m bp°ie "T, 2 , P oi. 



We have defined the reduced pole mass as m r = mi )PO i e m2,poie/ (^i, P oie + ^2, P oie) and x = 
1 — 4m r /(mi iPO i c + m 2i poic); when the two masses are equal, x — 0. The 's denote the 
coefficients of the QCD beta function: 

# ,> = 11 - f n,, /3^ = 102-|^, (5) 

£, and are given by 

t = ln(//r) + 7E , (6) 

a i ^' 

{ni) 4343 2 9tt 4 /1229 52Ca\ 100 2 ^ 

ai = h 367T + do G ni H nf- (8) 

2 lg S3 4 ^ 27 3 / ' 81 ' K J 

Up to C(a 4 m), the energy of a heavy quarkonium state X, identified by the quantum 
numbers n, I, s and j, is derived from the perturbative expansions of the energy eigenvalues^ 
of the above Hamiltonian: 

E x (n, ai ni) (/i) , m i)PO ie) = m 1)P<3 i e + m 2 ,p i e + KL,x a i Ul) it 1 ) > m i, P oie) , (9) 

^,^0*),^) = - ^^WVE^ 1 ^ Pk(Lnl), (10) 



1 The term H2, in its unitary-equivalent form, can be read from Ref. p~T| ] supplemented with the two- loop 
static potential calculated in Jl4| . 

2 The full formula up to 0{oqm) for the S state spectrum was derived in and later confirmed in [p~5| ; 
additional corrections necessary for the spectrum oH > 1 states can be found in [ fill and the formula for the 
unequal mass and / = case in fira. 
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where e = 1 is the parameter that will be used in order to properly organize the perturbative 
expansion in view of the C(Aq CD ) renormalon cancellation [J5J. P k (L n i) is a fc th -degree 

polynomial of L n \ = ln[3n/// (8a^ l \ii)m r ) }+Si(n+l) + -, and the harmonic sums are defined 

6 

q 1 

as = 2J 77^ ■ K is convenient to decompose the polynomials into renormalization-group 
invariant subsets: 

P = 1, (11) 
Pi = A^L^ + d, (12) 

3 Q(«i)2 r 2 , / l fl (r»i)2 . 1 n(ni) . 3 ( n ,) \ ( n; ) 



^ = 2& + (^/^ + t/T ; + ^o'ci) L nl + cT>. (13) 

ci and are given byj] 

ci = -4, (14) 
(no 16 tt 2 {2 S ( g + l)(l-x) + 3x} 87r 2 A(/, g ,j) x^.-C".)* , n 

C2 " TTn 6l0 + 9nl(l + l) (21 + 1) KM 

(ll+a;)7r 2 68 tt 2 473 9vr 2 33 Cs 9vr 4 /109 13 £ 



+ 77-7771 TT + + -7- + — Ti[ I "7777" + 



9n 2 9n (2/ + 1) 16 2 4 32 V 72 6 



(15) 



where X(l, s,j) represents the fine and hyperfine splittings for I > 0. If j = I ± 1 or m l polc 
m2, P oie (a; = 0), it is given by 



X(l, s,j) = -(1 - x)D s - (3 - x)X LS , (16) 



with 



n - /o (^- g ) 2 32 \ 2/(/ + l) g (s + l)-3X L5 -6X 2 5 
Ds =\ 3 —^- S ) = (2/-l)(2/ + 3) ' (17) 

X iS = (£ ■ S) = \ [j(j + 1) - 1(1 + 1) - s(s + 1)] . (18) 

If j = I and mi )PO ie ^ m,2, P oie (x > 0), the last term of Eq. (f|) induces a mixing between the 
s = and s = 1 states, and the total spin s is no longer a good quantum number (see, for 



instance, \T7\). In this case, the splitting is given 



X(l,±,j = 1) = 1± Jl + <ixl(l + l). (19) 



3 This corrects the formula given for c!j in Q in the situation when mi. po ic 7^ ??22,pole an d I > hold 
at the same time. Since, however, all the numerical analyses of @ correspond to the equal-mass case or to 
I = states, they remain valid. 

4 The corresponding energy eigenstate reads | ±) = [b \ s = 0) — A± | s — 1 ) ]/J \b\ 2 + Aj_, where b = 



2^J x I (I + 1) up to a convention-dependent phase. 
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The term v{n,l) is given by 5 



u(n, I) 



ir 2 1 n n + Z ! 

5 2 (n + Z) + — H — r 

8 2 2V ; 2 (ri- 



ll! S 



fc=i 



(n - Z + fc - 1)! 
(n + Z + fc)! fc 3 



| (n-Z-1)! "g 1 (2Z + fc)! (2fc + 2Z - n) 

k=l 



2(n + Z)! 



(fc-l)!(fc + Z-n) 3 



(20) 



It is understood that the last term is zero if n — Z < 2. 

Next we rewrite the series expansion of Ex in terms of the MS masses. This is done by 
expressing the pole masses mj )PO i e in terms of the renormalization-group-invariant MS masses 
Wi = ?rijMs( m iMs)- I n t ne theory with rii massless quarks only, we have up to 0(a 3 ): 



m 



i,pole 



4 ai ni) (mi) 2 (a^Hrrli) . 
i», < 1 + -e s V ; + e 2 s V ; < + e 6 

6 71 \ 7T / V 7T 



4 ni) 



(21) 



The coefficients 4 and d,2 are given in Appendix A . 

Note that the counting in e in Eq. (|I0D and Eq. ( [21] ) does not reflect the order in a s but 
the wanted renormalon cancellation [BJ. One way to understand this is to consider that in 
the sum of the pole quark masses and the static QCD potential, J2i m i,po\e + Vqcd(^), the 
renormalon cancellation takes place without reordering of power counting in a s || |4j] . The 
extra power of a s comes in the energy-level expansion when the dynamical variable r _1 is 
replaced by the dynamical scale (nlsj \ r' 1 | nlsj ) ~ Cpa B m r /n, where Cp = 4/3. Moreover, 
in order to realize the renormalon cancellation at each order of the expansion, it is necessary 
to expand m,j jPO i e and E^ n x in the same coupling |18|, |], [L9| ; we therefore express a 
in (pip in terms of a^ l \jj): 



rrii 



ft 



rrii) 



ft: 



(m) 



ft. 



rrii 



71 



ft'. 



K) 



i n 2 f M + ^ 



7T 



ln(A 



m,- 



• (22) 



Inserting Eqs. ( fffi) and (|2~1 ) into Eqs. ( |i~0|) and @, we get an expression for the energy 
levels of the heavy quarkonium states, which depends on /i, ft^ n ^(/i) and Mj, that we can 
organize as an expansion in e up to order e 3 : 



E x (n,ai ni \n),mi) 



_ , 4 a^'fe) 
^mJl + -e s 1 tJ +e 2 



4 n,) + e 2 



d 



(ni) 



71 



5 The infinite sum can easily be evaluated analytically in terms of £3, etc. for given values of n and I, e.g. 
by using Mathematica. 
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(23) 

« s (m ! )=Eq. (||) 



= m 1 + m 2 + E x l (1) (/i, a^K^^e + E% (2) (/i, c^(//), m,)e 2 

+^ (3) ( / i,a^( / i),m J )e 3 + ... (24) 

Since the counting in e explicitly realizes the order Aqcd renormalon cancellation, and since 
a s and rfii are short-range quantities, the obtained perturbative expansion ( ^4|) is expected 
to show a better convergence than Eq. fllPl). 

2.2 Physical parameters 

The input value for a s that we will use is 



= 

MS 


210^ 


MeV, 


= 

MS 


OQQ+30 


MeV, 


A®. = 

MS 


000+31 

000 -28 


MeV. 



af\M z ) = 0.1181 ± 0.0020. (25) 

Throughout this paper, we evolve the coupling and match it to the couplings of the theory 
with rti = 4 and 3 successively, by solving the renormalization-group equation perturbatively 
(analytically) at 4 loops (Eqs. (3) and (11) of Ref. |2T|).Q We obtain 

(26) 
(27) 
(28) 

In Ref. we calculated the bottomonium spectrum in the m c = case, i.e. directly 
from the above Eq. d24j).[] In particular, by fitting the theoretical values of the masses of 
the T(IS') and J/ip with the experimental ones, m& and m c were found to be 

m b = 4201+ig MeV (without charm- mass effects), (29) 
m c = 1237;^ MeV, (30) 

where we quote only the errors due to the uncertainty in a^ 5 '(Mz) given in Eq. Q2"5p. We refer 
to [U for an extended analysis of the other sources of uncertainty in the above determination. 
In [[|, with these values of the quark masses as input, other bottomonium masses have been 
calculated. The masses of the states 2 Si, 1P , 1P\, 1P 2 , 2P and 3 Si could be determined 
in a reliable way in the sense of (see also Sec. ||). 

6 We take the matching scales as mj, and m c , respectively. 

7 In Ref. Q major numerical results have been presented using the coupling ai (fJt) obtained by solving 
the 4-loop renormalization-group equation numerically. Here, we present the results computed using the 
above analytic coupling (corresponding to Table 2 column (i) of B). 
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Figure 1: Feynman graphs responsible for the leading charm-mass effects to the pole mass a) and to the 
Coulomb potential b). 

3 Charm-Mass Effects 

In this section we discuss how Eqs. @ and (]2"1~D get modified in the bottomonium case by 
including finite charm- mass effects up to order e 3 . These effects are generated by charm 
loop insertions and start at 0(af) in the first case [Fig. [!]&)] and at 0(a>D in the latter case 
[Fig. |l]a)]. Therefore, considering a charm with a non-zero mass modifies the right-hand side 
of Eq. @ (taken in the bottomonium case and with rti = 4 active and massless flavours) 
and of Eq. (|2T|) (taken in the bottom case and with n; = 4 active and massless flavours) 
respectively by an amount 

(6E bh ) mc = e\5E bl )£l + e 3 (8E bl )%l (31) 
(6m b ) mc = e\8m b )V+e%6m b yV. (32) 

The term (SE bb )^ c has been calculated for all quantum numbers in ||; (5m b )$ has been 
calculated in [|]; (5E bb )^ has been calculated in [|7| for the IS bottomonium level, but is 
unknown for higher quantum numbers; (Sm b ) $ has been calculated in ]7| in the limit of m c — > 
0, keeping only the (leading) linear contribution in the charm mass ("linear approximation"). 

In the following we will analyse these different contributions. Our main result will be the 
calculation of the dominant contribution to (5E bb )$ c (at the level of 95%) for all quantum 
numbers. Hence, we will be in a position to account for the finite charm-mass effects in the 
bottomonium spectrum up to order e 3 . 

3.1 Order-£ 2 effects 

The 0(e 2 ) effects in Eq. ( pip and Eq. ( |3"2"D stem from the graphs b) and a), respectively, of 
Fig. [T|. In [0| it was found that 



(Sm b 



(4)'-- 2 r 



,(i) _ m b I ar\ m b 



3 V 7T 



ln (£) + y - ( ln(£) + 2 K 



7 



0.015 



0.005 




1.2 1.4 1.6 l.S 



2 



. 02 



1.5 



2.5 3 



Fi gure 2: XTie /irsi figure shows (Smb)m\ (continuous line), (Snib)^ _>q (dashed line), {Smb)in\— »oo (dot- 
ted line), as a function of m c for rrib = 4.201 GeV. The second figure shows (Smb)ml (continuous line), 
(8mb)^ ^ (dashed line), (Smb)ml—>oo (dotted line), as a function of /i, when ai 4 \rn~b) is replaced byai 4 \ii), 
for rh~b — 4.201 GeV and m c = 1.237 GeV. The units are GeV. 



7T 



+ (1 + 0(1 + n Li 2 (-0 - - ln'tf) + HO Ml + + 7T 
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7T 



+ (1 - 0(1 - r) Li 2 (0 - - ln 2 (0 + ln(0 ln(l - - t 



(33) 



where £ = m c /m b . It is useful to define 



(8m b 
(8m b ) 



>m c ->a 
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(i) 

m c — >oo 



3 I 7T 



151 vr 2 13, (m c 

1 1 In — 

72 6 6 \m b 



+ \n 2 



m, 



(34) 
(35) 



where (5j7i&)mc-»o corresponds to the expansion of (8m b )$ for m c — > (up to 0(m;?) correc- 
tions) and (<fo7,&)mc-K>o corresponds to the expansion of (5mt)^ for m c — > oo (up to C(l/m 2 ) 
corrections). In the first plot of Fig. || we show (5mfe)^ as a function of fn c for m& = 4.201 
GeV (the value is taken from Eq. (|29|)). Here and in the following, if not specified differently, 
is calculated from af\M z ) = 0.1181 (central values of Eqs. (|26|)-(|2"8"|)). In the region 
of interest, 1 GeV < m c < 1.5 GeV, (8m b )^ turns out to be approximated by (8m b )^ c _^ 
reasonably well, while (Sirih)^^^ is far off. More specifically at the value m c = 1.237 GeV 
(£ = 0.294), taken from Eq. (30), we have 



(8m b )£l ~ 9.1 MeV, (5m&)g_>o ~ 10.5 MeV, (8m b ) ~ 18.7 MeV (36) 
The error of the "linear approximation" is about 15%, which agrees with the analysis done 



m 
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Since, in order to realize the renormalon cancellation, the coupling constant a^'(rfi b ) 
has to be expanded around the scale /i, in the second plot of Fig. |^ we show (5m fe )^, 
(Srrib)m c ^o an d {^ m b)m c ^oo as functions of /i at m c = 1.237 GeV and m b = 4.201 GeV, when 



ft 



(1) 



m b ) is substituted by oq 4 )(/j) in Eqs. (|33])-(|3"5|). Also this plot confirms that (5m b ) 



(i) 

m c — >0 



approximates (<5mb)W reasonably well, while (5m b 



is far off. 



-0.02 



-0 . 04 



-0 .06 



-0.08 




Figure 3: (SE)m\ (continuous line), (SE)!^ _, (dashed line), (SE)lnl^oo (dotted line), as a function of 
rh~ c formb — 4.201 GeV. Going down, the first set of lines corresponds to the IS state with /i = 2.446 GeV, 
the second one to the 2S state with fi = 1.065 GeV and the third one to the 35 1 state with /i = 0.724 GeV. 
Lines, which are not displayed, fall outside the plot range. The units are GeV. 

In IHl it was found that 



m b {C F af\p)) 2 af\p) f 3vr / . ^ . f ,\ . 

4n 2 3^ I Y n/? + T ( + } + (n + )(n ~ ~ T 

-7m Q(n + l)(2n + 1) + (n + Z)(n - Z - 1)) p 3 



+2 In 1-1 -2(^(n + Z + 1) - 



n-l-l 



(2n-l)! £ 



72n-l 



dp 



2n-l 



f n-l- 
i(2-P 2 



n + Z 
2Z + 1 + jfc 



p2(n-l-l-k) 



where p = 2nm t ./(m )) (7j7a^( 1 u)). Again, it is useful to define 



(5E t 



(i) 

bblmc^O 



-m. 



(1) 

m c —too 



m b {C F a^{p)fa^{p) f /2' 



(37) 



(3? 



4n 2 



3tt 



5 
3 



2^(n + Z + l)--0(l))},(39) 
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where (5E bb y r j l ' c ^ corresponds to the expansion of (5E bb )$ c for m c — ► (up to 0{m 2 c ) cor- 
rections) and {8Effi)$ corresponds to the expansion of (8E bb )$ c for m c — > oo (up to 
0(l/m^ +2 ) corrections). In Fig. [5] we show (SE bb )^ as a function of m c for m b = 4.201 
GeV. We plot the functions (from up to down) corresponding to the IS bottomonium level 
(ji = 2.446 GeV), the 2S bottomonium level (/x = 1.065 GeV), and the 35" bottomonium 
level (fj, = 0.724 GeV). The scales // correspond to the extremum of E bb (defined in Eq. 
(HP) at m b = 4.201 GeV. It is noteworthy that (8E bb )^ c _^Q does not depend explicitly on 
the quarkonium state (the dependence in Fig. |^ is only due to the different scales). In the 
region of interest 1 GeV < m c < 1.5 GeV, (8E bb )$ turns out to be very well approximated 
by (5^66)mc-^oo (f° r an the levels), while (5E bb )^ c ^ is far off. More specifically in Table p] we 
list for all the levels at the above scales the values of (SE bb )$ c , {5E bb )^ c ^o an d (<^&5)ml-).oo 
with m c = 1.237 GeV and m b = 4.201 GeV. The error of the "asymptotic approximation" 
_»cx>] * s extremely small. In the worst situation (the IS level), it is about 5%. In 
this case, however, the corrections are small (less than 1 MeV), so that they are beyond 
the accuracy of the present work. The "linear approximation" instead, at variance to what 
happens in the pole-mass expansion, is far off the exact values for all the states. This is in 
accordance with the general arguments of 0. Let us notice, as a general remark, that the 
"asymptotic approximation" for the energy levels, works better than the "linear approxima- 
tion" for the pole-mass expansion. We shall further comment on the above results in Sec. 
P. 



State 


n 




P 




(^Pf*)m^0 


( SE bb)ml^oo 


l 3 ^ 


2.446 


0.277 


1.59 


-0.0143 


-0.032 


-0.0136 


l 3 Po 


1.140 


0.428 


2.06 


-0.0210 


-0.076 


-0.0210 


1 3 P 1 


1.111 


0.437 


2.02 


-0.0221 


-0.079 


-0.0221 


l 3 P 2 


1.086 


0.445 


1.99 


-0.0232 


-0.082 


-0.0232 


2 3 S' 1 


1.065 


0.452 


1.96 


-0.0219 


-0.084 


-0.0211 


2 3 P 


0.726 


0.695 


1.91 


-0.0426 


-0.199 


-0.0424 


2 3 Px 


0.703 


0.733 


1.81 


-0.0490 


-0.222 


-0.0488 


2 3 P 2 


0.678 


0.782 


1.70 


-0.0581 


-0.252 


-0.0579 


S 3 ^ 


0.724 


0.698 


1.90 


-0.0405 


-0.201 


-0.0392 



Table 1: (5E bi )£l, (£%)m e -»o and ($ E bb)ml^oc for m b = 4.201 GeV and m c = 1.237 GeV; is 
calculated from A^L = 0.292 GeV at four loops. All dimensionful quantities are expressed in GeV. 

In Fig. £| we show (5E)$ (continuous line), and (8E)$ ^ (dotted line), as a function 
of \i for m b = 4.201 GeV and m c = 1.237 GeV for all the bottomonium states. Also this 
figure confirms that (SE)^^^ approximates (5E)$ very well for all the levels. 
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3.2 Order-e 3 effects 



The (Sirib)^ correction of Eq. (^) is not known exactly. However, following the calcula- 
tion of the two- loop mass effects in the static potential in [p2] , in [0 (^n^b)m c -*o has been 
calculated: 



-5(S + 21n2 ) + S (/ ' A + W }' (40 > 



where / 2 = 0.470 ± 0.005, 6 2 = 1.120 ± O.OlOfl j\ = (In A - ln& 2 )/(ln/ 2 - lnfo 2 ), h = 
(In A - ln/ 2 )/(ln& 2 - ln/ 2 ) and \nA= 161/228 + 13C 3 /19 - In 2. For the values m b = 4.201 
GeV and fn c = 1.237 GeV we have 

(5m 6 )2^ Q ~17MeV. (41) 

From the analysis of the previous section we may expect that also this value approximates 
(5mb)$ c with a relative uncertainty smaller than 20%. In [0] it is argued that the sum 
{Smb)m c ^o + {^ m b)m c ^o ma V indeed, approximate (5mb)^ c + (5mb)$ c with an uncertainty 
of 10%, while the precision of the sum (5rrib)ml + {^ m b)m c ^o is claimed to be worse. Finally, 
we notice that the series (^) shows no signals of convergence (compare with the figures of 
Eq. ©). 

The correction (5E b b)ml of Eq. ([H]) is known only for the IS* quarkonium state J7|. Using 
the analytic expression given in [0] at the values of the masses mj, = 4.201 GeV, m c = 1.237 
GeV and at the scale /1 = 2.446 GeV, we obtain^ 

(5E 1S )%] ~ -38.8 MeV, (^is)SUoc - "38-3 MeV. (42) 

This confirms that the error of the "asymptotic approximation" is, indeed, extremely small 
(in the case of Eq. (142]) it is about 1%). Considering that the IS" state is the state located 
furthest from the "asymptotic" limit, we may conclude that the uncertainty connected with 
the use of the "asymptotic approximation" in the energy expansion is for all levels below the 
accuracy of the present work and, hence, negligible. The explicit expression for (5E t 



bb 



,( 2 ) 



in terms of the pole mass is provided in Eq. (|68|) of |Appendix B| . 



3.3 Discussion 

In this section we interpret and discuss the above results. 

8 Somewhat larger errors are quoted in p3fl . In our analysis only the central values are used. 

9 We have used the MS masses in the expression given in M in terms of the pole masses. This introduces 
at order e 3 also a correction equal to 4a{ 4 \[i) / (3ir) (SE b i)m c (see also Eq. ( fl9| ) below). This correction 
amounts to ~ —1.7 MeV for the full (5E b i)ml and to ~ —1.6 MeV for the "asymptotic approximation" 
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a) The terms {8m^)^ and (5rrib)$ c are governed by the parameter £, which is the ratio 
between the charm and the bottom MS masses. This is natural, since these are the only 
scales involved in the pole mass expansion. Since £ ~ 0.3, which is smaller than 1, the small 
£, or m c , approximation is expected to work. On closer inspection of Eq. (|3"3"D, we find 

(6m b )% - (Jm»)£U = _6_ + Q{e) ^ _ lg%j (43) 
(Sm b )m c 71 

which practically accounts for the entire difference between the first two terms in (|36|) . 
Therefore, we may say that, concerning the pole mass expansion, the charm mass can be 
considered small and we are close to the situation of four active and massless flavours. 

b) The terms {5E bb )^ c and (5E bb )$ c are governed by the parameter p. This is twice the 
ratio of the charm MS mass and the typical momentum of b and b, which is the relevant 
scale of the bound state. This number turns out to be typically larger than 1 (see Table |l|). 
Therefore, an expansion in small p, or m c , is out of the question and only an expansion for 
large p, or m c may work. From Eq. ([37]) we have 



(5E bi )£l - g£^)gUc 1 1 (u) 

(SE bb )g\ ~ p«+ a 2V>(n + Z + 1) - 2^(1) + 5/3' 1 ] 

It follows that, even if the expansion parameter p is not particularly large, the "asymptotic 
approximation" turns out to work very well (for all the states), since the leading discarded 
term is suppressed as l/p 2l+2 , i.e. by at least two powers in 1/p. The numerical factors 
contribute for extra suppressions. This explains why the data of Table [I] are reproduced so 
well by the "asymptotic approximation". Notice that, i) states with higher I are expected 
to be reproduced better by the asymptotic formula than states with lower I (this can be 
checked explicitly in Table |I] by comparing the data for the / = states with those for the 
I = 1 states, and in Fig. £|); ii) states with higher n are expected to be reproduced better 
than states with lower n (as long as we are in the perturbative regime), since the parameter p 
grows like n/a^\p) (even if the growing of a^'(p) partially compensates the growing of n). 
In this sense the conclusions drawn in Sec. |3.2| on the relative size of (SE bb )^ c with respect 
to (SEjji)^^^ for the IS state are expected to hold even better for the other (higher) states. 

We may conclude that, where the energy-level expansion (in terms of the pole mass) is 
concerned, the charm mass can be considered large with respect to the typical scale of the 
bound state, so that the charm effectively decouples. Hence, charm-mass effects can be taken 
into account in the energy-level expansion of the bottomonium system (in terms of the pole 
mass) in a very effective way by considering it in the situation with three active and massless 
flavours onlyQ. We have 

(SE bb ) mc 9 f % (tf-EtfW-oo = El m bb (p, ai 3 \p),m btPole ) - E^ bb (p, c^ 4) (/i), m 6iPole ), (45) 



10 It is worth while to stress that the situation here is the exact opposite with respect to that of the 
pole-mass expansion in terms of the MS mass discussed in a). 
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where ,7 are the bottomonium binding energies in terms of the bottom pole mass (Eq. 
(|ITi|)) with rii active and massless flavours only. Equation ( f4"5] ) is the key result of our analysis. 
It is expected to be valid at any order and it states that the charm decouples in the energy- 
level expansion (in terms of the pole mass) of the bottomonium system. Explicit expressions 



for Eq. (|45D up to order e are given in [Appendix B 



3.4 Cross-check 

We have made a non-trivial cross-check of the major result, Eq. (|45|) . We compared the total 
static energy of a bb system, E tot {r) = 2m bpole + Vqcd(^), up to (D(a 3 ), in the case where 
the full charm-mass corrections are included in VqcdM; the static potential in perturbative 
QCD, and in that where Vqcd(^) is evaluated in the limit m c — > 00. To be more precise, we re- 



analysed the total energy E tot (r), studied in two hypothetical cases in [p4}| . We have included 
the full charm-mass corrections to the static QCD potential using the formulas obtained in 
P2| , 0. Then we compared the total energy with the one calculated using Vqcd(^) m the 
limit m c — > 00, i.e. with rii = 3. [In both cases, (5m b )^ c is approximated by (5m b )^ c ^ .] 
The total energies in the two cases agree with each other to a good accuracy: the difference is 
3 MeV at r = 0.5 GeV" 1 and of order 0.1 MeV in the range 1 GeV" 1 < r < 3 GeV" 1 . Since 
Etot{ r ) determines the bulk of the bottomonium spectrum ||, and since all the charm-mass 
effects are included in it, this confirms the validity of our approximation. Note that the 
cross check was made after cancellation of the leading renormalon contributions (expressing 
m^poie in terms of which adds a non-trivial point to the analyses done in the previous 
sections. 



3.5 Level expansions in terms of MS masses 



A way to implement the above results in the bottomonium level expansion in terms of the 
MS masses is to modify Eq. (|23|) into 



E 



bi, 



2fn b <l + -e 



A^af\m b ) t A (af\m b ) ^ (4) ( ^ f a< 4 > (m b ) \ 3 rf(4 ) 



+ e 



7T 



7T 



TV 



+£ 2 2(* mi )S + £ ^2(^)g)-i( c '' tt - 4,( '' ) ' 



n 



(8m b ) 



(i) 



a s (m h ) =Eq. (22) 
ai 4) =Eq. © 



(46) 



2m b + E$ (//, a® (fi) , m c , m b )e + E$ (//, a^ 3) (ji) , m c , m b )e 2 
+E%\n,ai 3 \ti),m c ,m b )e 3 + ..., 



(47) 
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where a^'(rn b ) in Eq. ( fj6|) is understood as expanded everywhere around the scale p ac- 
cording to Eq. and around the coupling with 3 active massless flavours according to 
the following equation [21]: 



«i 4) G") = a^(p)ll+s 



(3), 



a^(p) 



3tt 



2 r 



7T 



9 



In 2 



12 Vm r 



11 
72 



(4* 



Both expansions are needed in order to have the same coupling in the pole mass and in the 
binding energy expansion so that the (9(Aqcd) renormalon cancellation is made explicit. 
Finite charm-mass corrections affect the pole-mass expansion and are explicitly written in 
the last line of Eq. ©. 

Another possibility is to express also the binding energy at four massless flavours and to 
calculate the corrections to it. In this case Eq. (|23| ) is modified into 



E, 



i>b 



2l7lh 



1 , 4 c aW(m fe ) | g2 ( af\m b ) 

3 7T V 7T 



4 4) + e 



(I 



71 



('1) 
2 



b'm,bb 



p, af\p),m b I 1 + 



4 af\m b ) 



+ e 



2 / af\m b 



71 



71 



(4) 



+e^2(5m b )W+(8E bl )mm b ,p) 



Mr™, 



+e 3 ^ 2(5m fe ) 



(2) _ i 



1 /CVo^Y 



4 V 



n 



(5m b ) 



(i) 



H6E bb )W(m h p) + ~?^(5E bk j£l(m h p) 



J Qs (m 6 )=Eq. (22) 



= 2m b + Ej£ (//, a^ 4) (//) , m c , m b )e + ^Jf (/i, c^ 4) (p) , m c , m 6 )e 2 



v 6& 

+ E lf 0"> a s 4) ()") > m c , m b )e 6 + . . . , 

where we have made use of the fact that p = p(l + O(a^)), and where a^\fn b ) in Eq. ( |49"D 
is understood as expanded everywhere around the scale p according to Eq. (p2[). The first 
two lines of Eq. ( f4Tj|) exactly correspond to Eq. (|23|), when calculated in the bottomonium 
case with four massless flavours. 

Let us summarize our present knowledge of the different pieces of Eqs. (f|6]) and (|49|) . 

A ) E^ n bb - are the bottomonium binding energies in terms of the MS masses with ni active 
and massless flavours only (i.e. m c = 0). They can be read from Eq. (|T0|) . A study of these 
contributions alone can be found in Ref. @. 

B) (Sm b )$ c is given in Eq. d33|). If we expand a^(jn b ) around the scale p, we then have 
to apply the replacement of Eq. fl2"2"|). If we expand a^'(p) around a^(p), we then have to 
apply also the replacement of Eq. (f48[). The quantity (8m b )^l is known only in the small m c 



(49) 



(50) 
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The use of (5mt,)m-+o instead of (8mb)$ should account for 



it is claimed that the simultaneous use of (Jmj)^'^ and (5m b )^^ 



limit and is given in Eq. (|4 
80% of the effect. In § 
should account for 90% of the actual effect. 

C) (5E bb )^l(m b , p) is given in Eq. (|37|). The large m c limit of this expression is given 
in Eq. (|39| ) and practically accounts for the full effect. The term {SE^)^ (m b , p) has been 
calculated here in the large m c limit and is given in Eq. (^). For the reasons discussed in 
Sees. ^3] and we believe that Eq. (|68D accounts for practically the full effect. 

Since Eq. (f46|) accounts explicitly for the decoupling of the charm quark at all orders in 
the binding-energy expansion (in terms of the pole mass), we expect results in this case to 
be more stable and reliable than by using Eq. (|49|). This will be confirmed by our numerical 
analysis. Therefore, we shall use Eq. fl46|) as our reference calculation (Sec. f4.2|) . Equation 
(||9|) will be used in our analysis of theoretical uncertainties (Sec. fOl) . Moreover, in our 
reference calculation of Sec. JO] we will approximate both (5m b )^ c and (5m&)^ c by means 



of the "linear approximations" (5mb)m c ^o an d (^ m fe)mc->o- This proves to be slightly more 
stable than the use of the full correction (Smb)^- We take into account the full correction 
in the analysis of the uncertainties in Sec. 



v(2) 



We conclude this section by giving the expression, which substitutes Eq. (£3|) in the B c 
case, when finite charm-mass effects are taken into account: 



E 



Br: 



m c < 1 H — e 



4 ai 3 Hm, 




+E'L A I /»,«?"(/<) 



n 



v m c + m b J 



ai i] =Eq. (P 



= m c + m b + Egl (ji, a? } (p) } m C} m b )e + Eg] (//, a? 5 (//) , m c , m b )e 2 
+EP c (p, ai 3) (p),m c ,m b )e 3 + ..., 



a s (m c , b ) =Eq. (22) 
(51) 

(52) 



where a4 3 ) (m c ) and (m b ) in Eq. (|5T|) are understood as expanded everywhere around the 
scale /i according to Eq. (|22"D and a s "(/x) around a^ 3 )(/i) according to Eq. fl48|). 
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4 Numerical Analyses 

In this section we examine numerically the series expansions of the bottomonium spectrum, 
Eqs. (|D and flU), and of the B c mass, Eq. (P). 



4.1 Scale-fixing procedures 

The quarkonium mass Ex depends on the scale fi, because of our incomplete knowledge of 
the perturbative series. We will fix the scale fi for each state X in two different ways. 

A) We fix the scale \i = /i x by demanding stability of Ex against variation of the scale: 

= 0. (53) 



E x (fi,ai ni \fi),mi) 



B) We fix the scale jj, = [i x on the minimum of | EjP | : 



d_ 



[Ef{^af l \^i)? 



0. 



(54) 



When we do this, we expect that the convergence properties of the series become optimal, 
and that the scale becomes close to the inverse of the physical size, ax, of the bound state 
X (defined as in Ref. 0). If the scales fixed by Eq. (|53| ) and Eq. ( |54]) evidently do not 
fulfil these expectations, the theoretical predictions obtained in this way will be considered 
as unreliable. This typically happens when the coupling constant becomes larger than 1. In 
the next section we will check if \E%\ > \Ex \ > \Ex\ &t the chosen scales. Notice that 
the criterion B) of Eq. fl5"4|) practically always fulfils this request since min{|£' x 3 ^|} = in 
all the considered cases. This means that the results obtained using criterion B) have to 
be taken with some caution. For this reason, in the following we will use them, rather than 
as an independent way to calculate the heavy-quarkonium spectra, as an independent check 
of the results obtained using criterion A) as well as a means to determine the uncertainties 
related to the fixing of the scale fix- Comparisons of the chosen scales and the bound-state 
sizes have been given in Sec. 5 of and we do not repeat them here. 



4.2 Numerical results 

Here, we present results from both scale-fixing procedures on Eq. (|46| ) for the bottomonium 
spectrum and on Eq. ( [5T] ) for the B c mass. We calculate both (Smb)$ c and (Srrib)^ by 
means of the "linear approximations" {Smb)m c ^o an d (^ijm^o- The coupling constant a s 
is defined as in Sec. [T21 The uncertainties on af\M z ) [or on Aj^g, see Eqs. (p6|)-(p8|)1 are 
the only ones that we take into account here. Other sources of uncertainties and in particular 



those related to the finite charm-mass effects will be discussed in Sec. 4.3 



17 



btate A 


r-iexp 




771CXP 771 


pi 






Vx 


„ (3) 1 ,,A \ 


I (1 bt) 


y.4oU 


O A RC\ 

y.4bU 


U 


U.000 


U.zUo 


U.UUb 


Q 1 A 

z. 14 


U.zoo 




y.obU 


n nntx+75 

9.995_ 62 


-0.135 +62 


1 £ O A 

l.oo4 


U.1U1 


—0.021 


1 no 
l.Uo 


n a 

o.4oy 


I (1 fx) 


n coo 
y.oyo 


10.004_ 63 


n 1 11 —78 
-0.111+63 


1.004 


U.Uol 


— U.UZZ 


l.Uo 


U.4oo 


Y/i3d \ 
1 (PI2) 


9.913 


10.012165 


r\ AAA- 81 

-0.099+65 


1.591 


0.063 


—0.022 


1 C\0 A 

1.034 


0.477 


T fo3 C \ 

I (2 bi) 


lU.Uzo 


1U.Uo4_ 75 


-0.061+ 7 5 


1 fii s 
l.olo 


U.Uyb 


— U.U1U 


1 no 
l.UZ 


U.4oo 


Tf9 3 P ^ 


1 n oqo 

HJ.ZoZ 


in y o+196 


n « — 196 


9 /I91 


— u.000 


n 1 A9 
U. 1UZ 


A 77S 
U. ( ( 


n 71 n 
U. i 1U 


T(2 3 Pi) 


10.255 


10.5641?™ 


-0.309;?™ 


2.472 


-0.404 


0.116 


0.770 


0.726 


T(2 3 P 2 ) 


10.268 


10.578+?! 


-0.310+11 


2.518 


-0.449 


0.129 


0.762 


0.740 


T(3 3 5!) 


10.355 


10.645i? 6 | 


-0.290;?68 


2.472 


-0.348 


0.140 


0.770 


0.726 


T(4 3 5x) 


10.580 


* 


* 


* 


* 


* 


* 


* 


^(1%) 


6.4 ±0.4 


6.307±| 


0.1±0.4 


0.675 


0.188 


0.017 


1.62 


0.334 



Table 2: Theoretical predictions for the bottomonium and B c masses from Eqs. ^h\ ) and (j5Jp, respectively 
((8mb)m'a^ — (^ m b)m'^*o)i when the scales \i x are fixed through Eq. (|5^). The quantities ar ^ the order 
£■? contributions to the spectrum and are defined by Eqs. and \5^). The c- and b-quark MS masses 

are fixed on the experimental values of the J/tp and Y(IS) masses, respectively. The coupling constant 
ai 3 ^ has been defined according to Sec. The uncertainties in the third and fourth columns refer to the 

uncertainties in ai 5 \Mz) only. All the other data refer to af'fM^) = 0.1181 and to the MS quark masses 
fixed on the central values of Eqs. and (I5T). All dimensionful numbers are in GeV. 



Scale-fixing procedure A ) 

Since we expect the ground states of the bb and cc systems to be the states least affected by 
non-perturbative corrections, we fix and m c through the two conditions 

E r{1S) (lJ,i,ai 3 \fi x ),m c ,m b ) = E c ^ s) = 9.460 GeV, (55) 
£ w (/4,af (/4),m c ) = = 3.097 GeV, (56) 



where the experimental values of the vector ground states have been taken from [20]. We 
assume, for the moment, that this identification is not affected by non-perturbative correc- 
tions. From Eqs. fl55|), fl5"6|) and fl5H|) we determine [i x (see Table |2|), and the b and c MS 

masses: 

m b = 4190;?^ MeV, (57) 
fn c = 1237+^j MeV. (58) 

These values are in good agreement with the estimates based on T |25|, [y| [7j and 
charmonium |27| sum rules respectively. The charm mass given in Eq. ([58]) corresponds to 



the one quoted in Table 2, column (i) of Ref. (reported also in Eq. (|30"D of the present 
paper), since it is not affected in our analysis by massive quark- loop effects. The bottom 
mass given in Eq. fl5"7|) is new. The errors refer to the uncertainties in af\Mz) only. 
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Using these masses as input, we calculate the energy levels of other observed bottomonium 
and B c states. In Figs. [5|-[7] we display the \i dependence of the different energy levels E x 
measured from IrfTb and nib + m c , respectively (at different orders in e) before scale fixing. 
The levels, after scale fixing by Eq. (^), are given in Table |2|. Theoretical predictions that 



we consider unreliable, in the sense specified in Sec. |4.1| , are not displayed and indicated 
with a star (*). Generally, for states that we consider reliably calculable in the perturbative 
approach, the scale dependence decreases as we include more terms of the perturbative 
series. For states whose predictions we consider unreliable, the series would become much 
more convergent if we chose a scale different from (typically larger than) fix ■ The theoretical 
prediction for B c is consistent with the experimental value, although the experimental error 
is large. It is also in agreement with the determination of [Rfl. 
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Figure 5: EQ (dotted line), E^ + E^ (dashed line) and E^ + + E^ (continuous line), as defined 
in Eq. fcfy ) ((Smb)m^ — (Srnb)^' 2 ^), versus fi for the states lSi, 1-Po; an d XP2- The parameters 
are = 4.190 GeV, m c — 1.237 GeV and A^l. = 0.333 GeV; ji x obtained from the minimal sensitivity 
prescription, Eq. ft53), is explicitly shown. The units are GeV. 
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in Eq. (|5qj versus \i for the B c . Other notations are the same as in Fig. 0. 
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Table 3: Theoretical predictions for the charmonium, bottomonium and B c masses, when the scales /z^ are 
fixed through Eq. The other conventions are the same as in Table [|. 
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Scale-fixing procedure B) 

By fixing the bottom and charm MS masses on the experimental values of the T(1S) and 
J /if) masses via the scale-fixing procedure of Eq. (0), respectively, we obtain 

m b = 4190+^ MeV, (59) 
m c = 1241+JI MeV. (60) 

These values are in agreement with the estimates based on the procedure A), given in Eqs. 
(|57|) and ([58]). The errors refer to the uncertainties in a$>(Mz) only. Using these masses as 
input and fixing the scales through Eq. ( E4j ) , we calculate the energy levels of other observed 
bottomonium and B c states. The levels are given in Table |||. Reliable theoretical predictions 
agree well with those given in Table ||]. 

Let us compare the present results for the bottomonium states with the corresponding 
results obtained in the m c = case in Table 2 column (i) of 0]. The value of the bottom 
MS mass m& determined from the T(15) mass is lowered by about 11 MeV by the inclusion 
of the non-zero charm-mass effects, for the same input value of af\Mz). The masses of 
the n = 2 bottomonium states are raised by about 70-100 MeV by the charm-mass effects; 
the masses of the n = 3 states are raised by about 240-280 MeV by these effects. These 
tendencies can be understood as originating from an increase of the binding energies in these 
systems: if \i x ,B < m c , at the corresponding scale becomes larger than a^; hence, the 
higher the energy level is, the more its mass is increased. For this reason, the value of 
the input af\M z ), which reproduces the whole level structure of the experimental data 
optimally, becomes smaller when we include the charm-mass effects: in the m c = case, 
ctf\Mz) = 0.1181 reproduces the experimental data fairly well, if we take an average of 
the S'-wave and P-wave levels corresponding to each principal quantum number n, while 
the agreement is better for (Mg) = 0.1161 when the charm-mass effects are included. It 
also remains that the predictions for the S—P splittings and the fine splittings are smaller 
than the experimental values, as noted in p| . However, it should also be noted that for these 
observables the perturbative expansions considered here include only tree-level contributions. 
As for the B C (1S) mass, it is reduced by about 10-20 MeV by the charm-mass effects, when 
m&, m c and a s (Mz) are taken as input. 

4.3 Error estimates 

In this section we take the results listed in Table Q as our reference calculation and discuss 
different kinds of uncertainties that may affect them. Besides non-perturbative corrections, 
there are four kinds of uncertainties, which we list below. 

1) Uncertainties related to the charm-mass effects discussed in this paper. Since the 
charm-mass effects are the main subject of the present paper, their errors will be 
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Uncertainties of charm-mass effects 
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Table 4: Differences between the determination (i) (ii) and (iii) discussed in the item 1) of Sec. ^.i and 
the theoretical predictions of Table [|. The central values correspond to ai 5 \Mz) = 0.1181, the upper values 
to ai 5 \Mz) = 0.1201, the lower values to OLs^(Mz) = 0.1161. The scale is fixed according to the procedure 



A) of Eq. ([530 . Those values corresponding to unreliable theoretical predictions are indicated with a star (*). 
All dimensionful numbers are in MeV. 



discussed in some detail. In evaluating the effects of charm-mass loops in the pertur- 
bative expansion of the quarkonium levels, we have considered different options (see 
Sec. H). Here we estimate the differences between these options and the reference 
calculation in the previous section. This also shows that the reference calculation is 
the most reliable in the sense discussed at the beginning of Sec. £|. In Table f| the 
following error estimates are listed (the scales \i are fixed by Eq. (^3|); the quoted val- 
ues correspond to af (Mz) = 0.1181 ± 0.0020): (i) The difference between the results 
obtained from Eq. ([46]) (Eq. ([Hp for B c ) when the full expression of (5rrib)$ , as given 
in Eq. fl33|), is used instead of the linear approximation (5m b )^ c ^ and the results of 
Table (ii) The difference between the results obtained from Eq. fl4"9p when (Srrib)^ c 
is approximated by ((Jmj)^^ and the results for the bottomonium system of Table ^. 
(iii) The difference between the results obtained from Eq. ( |49| ) when the full expression 
of (SmbYml is used and the results for the bottomonium system of Table 0. 

2) Uncertainty of a^\Mz). This has been taken into account in the results presented in 
the previous section and listed in Tables |] and [3|. 

3) Uncertainties from higher-order corrections. These have already been considered in 
and may be estimated as follows (see Table [5]): (i) The difference between the theoreti- 
cal predictions computed using a s (/z) as obtained by solving the renormalization-group 
equation numerically at 4 loops and the definition of a s (/i) used in the previous sec- 
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Uncertainties from higher-order corrections 
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Table 5: Variations of the theoretical predictions of Table when the uncertainties 3) (i) and (ii) of Sec. 



4-S are separately taken into account. The number in bracket shows a naive estimate of uncertainties from 
different e-expansion schemes. Those values corresponding to unreliable theoretical predictions are indicated 
with a star (*). The input parameter is a s {Mz) = 0.1181. The scale is fixed according to the procedure A) 
of Eq. All dimensionful numbers are in MeV. 

tions and specified in Sec. |2.2| . (ii) The contributions ±|P^| from Table || (±|P^|/2 
for Srrib). We have not explicitly considered here uncertainties related to different im- 
plementations of the e-expansion scheme. Indeed, one could imagine different ways 
of counting in e terms that are not directly related to the renormalon cancellation in 
the pole mass and in the static potential, such as relativistic corrections or spin-orbit 
and spin-spin interaction terms. A naive estimate, obtained by subtracting from the 

/o\ _____ 

definition of E x given in Eq. (H7]) all the terms mentioned above, gives in the m_ case 

/o\ 

± | P^ (after reshuffling) | /2 = ±25 MeV. However, further investigations are needed 
on this point. In || also the difference between the theoretical predictions computed 
using the 3-loop and the 4-loop running coupling constants has been considered. The 
size of these uncertainties is typically smaller than those listed here, and, therefore, 
they are not considered. 

4) Uncertainties in the scale-fixing procedure. These are also part of the higher-order 
corrections and may be estimated as follows (see Table [|): (i) The difference between 
the theoretical estimates of Table |3] obtained by fixing fix = f^x v i& the condition ( |5"4"D 
and the results of Table 0, obtained by fixing fi X = f^x v * a the condition (^). (ii) 
For comparison with the error estimates given in conventional analyses, we consider 
the (maximal) variations of m_ and Ex when we fix the scale as ±10% of the minimal 
sensitivity scale: /i = jj, x x (1±0.1), where \i x is taken from Table [| This is numerically 
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Uncertainties in scale fixing 
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Table 6: Variations of the theoretical predictions of Table when the uncertainties 4) (i) an d (H) of Sec. 



4-o are separately taken into account. Those values corresponding to unreliable theoretical predictions are 



indicated with a star (*). The input parameter is a s (Mz) — 0.1181. All dimensionful numbers are in MeV. 



close to the scale- fixing uncertainty evaluation done in 12 . 



A summary of the above uncertainties is given in Table [7[ The figures are obtained as ± the 
absolute value of the maximum of the corresponding uncertainties. In the column labeled 
"h.o. corrections" the maximum is taken from the uncertainties of Tables |5| and []. Note 
that in Table |5| we have not calculated the error associated with different implementations 
of the e-expansion for levels higher than the ground state. 

4.4 Charm-mass corrections 

In Table |]we display the pure charm-mass corrections of order e 2 , 

2(5m b £Lo + (61) 

and of order e 3 , 



4 \ n j 3 7r 

(62) 

of Eq. ([E]) , where the energy is expressed in terms of the 4-flavour coupling. The full energy 
levels can be easily obtained from the data displayed in column (ii) of Table |] and from 
Table |[ Table |8| shows that also the disentangled finite charm-mass corrections exhibit a 
reasonable convergence, in correspondence to the reliable predictions of Table |2|. 
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Table 7: Summary of the uncertainties discussed in Sec. \4 $ - The figures are obtained as ± the absolute 
value of the maximum of the corresponding uncertainties. In the column labeled "h.o. corrections" the 
maximum is taken from the uncertainties of Tables [| and All dimensionful numbers are in MeV. 

5 Conclusions and Discussion 

We have computed the bottomonium spectrum within the framework of perturbative QCD, 
including the non-zero charm-mass effects. The effects of a charm mass m c > as compared 
to the m c = case are to increase the level spacings; the effects are larger among the higher 
levels. They can be understood as follows: the effective coupling becomes larger at the 
relevant scale when the decoupling of the charm quark is incorporated; consequently, the 
binding energy increases. Since we fixed m& = m^ IS (m^ IS ) on T(IS') in our analysis, the net 
effects are to decrease % by about 11 MeV and to increase the n = 2 and n = 3 levels by 
about 70-100 MeV and 240-280 MeV, respectively. 

We have analysed the uncertainties of our predictions within perturbative QCD. The 
uncertainties originating from the error of the input af\Mz) [Sec. [p: 2)] are as large as (in 
most cases larger than) other uncertainties, reflecting a high sensitivity of the bottomonium 
spectrum to (M^). The uncertainties from unknown higher-order corrections [Sec. |4.3| : 3) 
and 4)] are consistent with the estimates based on the next-to-leading order renormalons, i.e. 
are numerically of the same size as Aqcd x (o-xAqcd) 2 : if we approximate I /ax — Hxi take 
the values of Table ||, and Aqcd = 300-500 MeV, we obtain for the IS state a contribution of 
order ±(5-30) MeV, for the n = 2 states a contribution of order ±(20-130) MeV and for the 
n = 3 states a contribution of order ±(40-220) MeV. The uncertainties of the charm-mass 
effects have been discussed in Sec. fO|: 1). 

If we take a lower value (~ 0.1161) of the input af'(Mz) within its present error, the the- 
oretical predictions for the bottomonium spectrum are in agreement with the experimental 
data within the next-to-leading renormalon uncertainties. Considering the high sensitivity 
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Table 8: Pure charm-mass correction of order e 2 and e 3 to Eq. (last three lines). The input parameter 
is a s (Mz) = 0.1181. All quantities are in GeV. 

to a^(M z ), the agreement (inside the present world-average values) seems to be quite non- 
trivial. Thus, the data at our disposal would confirm the analysis done in 0: (1) the bulk 
of the bottomonium spectrum is accessible by perturbative QCD up to the n = 3 states 
(at least in the scheme of Table ||D; (2) non-perturbative contributions may be of the type 
associated with the next-to-leading renormalon as estimated above. 

If the true value of a^(M z ) turns out to be close to its upper value of 0.1201, however, 
independent non-perturbative effects, which cannot be absorbed into the next-to-leading 
renormalons, should exist, with magnitudes of at least about 60 MeV for the n = 2 states 
and at least about 260 MeV for the n = 3 states. In this case, the non-perturbative effects 
should decrease the higher levels, resulting in a tendency opposite to that of the effects 
usually considered by local gluon condensates or by a linear confining potential. Unless 
a^(Mz) is known more precisely, we cannot constrain the non-perturbative effects more 
stringently with our method. On the other hand, once (M^) is measured more precisely 
in the future, it is certain that we will be able to investigate the non-perturbative effects on 
the bottomonium states very much in detail, since there are a number of energy levels that 
can be computed reliably in the present perturbative method. 

Yet another possible approach is to take advantage of the high sensitivity of the bot- 
tomonium spectrum to af (M^) and to determine both af\Mz) and mb from a fit to the 
experimental data, assuming that non-perturbative effects can be absorbed into the next- 
to-leading renormalon uncertainties. Then we may compare the value of a^(Mz) with 
those obtained from other measurements, for a consistency check of the assumption. It is a 
meaningful method to study the nature of the non-perturbative effects. 

Therefore, according to our observation that the perturbative predictions, within uncer- 
tainties, are consistent with the experimental data, we present the value of the bottom quark 
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MS mass as determined from the Y(IS') mass in our approach 

m b = 4190 ± 20 ± 25 ± 3 MeV, (63) 

where the first error refers to the uncertainty in a^(Mz), the second error to the higher-order 
corrections, and the third error to the uncertainty of the charm-mass corrections. It is useful 



to compare the above result with a recent estimate done in p8fl , mj = 4210 ± 25 ± 90, where 



the renormalon cancellation has been implemented in a scheme different from the e scheme 
adopted here. Taking into account that in p8| charm- mass effects have not been included 
in the calculation (only added to the errors) the central value of 4210 MeV is consistent, 
within errors, with our central value of 4190 MeV [or, with the value of Eq. The 
error ±25 MeV refers to uncertainties in af'(Mz) and corresponds to our ±20 MeV. The 
error ±90 MeV refers to charm-mass effects, higher-order perturbative and non-perturbative 
corrections. It corresponds to our errors ±25 ± 3 MeV. Subtracting from ±90 MeV a ±(10- 
20) MeV estimate of the charm-mass effects that have not been calculated explicitly in ||28|| , 
but have been included in our analysis, it remains a ±(40-50) MeV difference between the 
two error estimates. The difference is given by the size of the non-perturbative corrections. 



These are guessed to be ±50 MeV in p8[ . On the other hand, in the present work, the 
observation that the (perturbative) determinations of the n = 2 and n = 3 bottomonium 
levels agree within the present uncertainties with the experimental data suggests this to be 
also the case for the n = 1 level, i.e. for nib- In other words, non-perturbative effects would 
be absorbed into the present uncertainties. Hence, higher bottomonium levels provide the 
additional physical input needed in order to constrain the non-perturbative corrections on 
the b mass. We also mention that in the analysis of 0, where m& has been extracted from 
the T sum rules including charm-mass effects, the author obtains = 4.17 ± 0.05 GeV. 
This figure is consistent with our estimate. 

Let us, finally, remark some general tendencies of the perturbative predictions. First, 
evolving a^ ni \/j,) from a^(Mz) using the 3-loop running instead of the 4-loop running, is 
theoretically still consistent within our present accuracy. In this case, the level spacings of 
the spectrum become smaller, since the 4-loop coefficient of the beta function is negative. 
Secondly, if we compare the analytic running-coupling constant (defined in Sec. [2.2| and 
used throughout this paper) and the numerical running coupling constant (used in 0), the 
former is larger at relevant scales. Hence, the level spacings become smaller if we use the 
latter coupling. These variations of the predictions are certainly higher-order uncertainties. 
Nevertheless, it is surely worthwhile to further examine the quality of the theoretical predic- 
tions carefully and to seek for most reliable theoretical predictions at the current reachable 
accuracy. 
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Appendix A Pole mass expansion coefficients 



Here, we give the analytic expressions for the coefficients 4 and 4"^ that appear in Eq. 
([2l"l). The coefficient 4 has been calculated in ||, while the analytic expression of 4"^ 
can be derived from the result of [ 30|] , 11 the renormalization-group evolutions of a^'^/x) and 



m Ms(/ i )' anc ^ the matching condition fTilfn 



{ni) 307 tt 2 7r 2 ln2 Cs^ '71 ^ 



4 ni) 



32 3 9 6 V 144 18 / 

13.4434- 1.04137 nj, (64) 

8462917 652841 vr 2 695 vr 4 575 tt 2 In 2 



93312 38880 7776 162 

22 7r 2 ln 2 2 55 In 4 2 220Li 4 (±) 58 Ca 1439vr 2 C 3 1975 C 5 



81 162 27 27 432 216 

/ 231847 991 vr 2 61 tt 4 Il7r 2 ln2 2vr 2 ln 2 2 In 4 2 

+ni \ 23328 ~ 648 + 1944 81 + 81 + ~8T 

8Li 4 (|) 241 ( 3 \ 2 / 2353 13tt 2 7 ( 3 \ 

27 72 J Ul \ 23328 324 54 ) 

~ 190.391 - 26.6551 + 0.652691 n]. (65) 

Appendix B Finite charm-mass effects in the 

energy expansion 

In this appendix we give the expression of (5E bb )!^ ^^(rrif,, p) for a bb state of quantum 
numbers n, I. From Eq. fl4"5|) we get after explicit calculation (we write here explicitly the 
mass arguments, which are the 6-pole mass and p = 2nm c ^ \ e /(m b ^ \ e CFCt { f ) (p)))'- 

(££&&)m c -*oo = e 2 (8E bb )^ 00 (m b , po i e , p) + e 3 (8E bb )^ 00 (m btP0 i e , p), (66) 
<«fc>£U>*-../>> = ra, r M\,))^y 

*{i h G)-M(*» + , + 1 > + , »)}- (67) 



11 This constant was obtained numerically in |29fl in a certain approximation. 

12 This relation coincides with Eq. (14) of |50|, which is given numerically (indirectly through /3o and j3\). 
Note that, in the other formulas of |30|, the coupling of the full theory is used. 
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Equation fl5"T| ) corresponds to Eq. 
limit with Eq. (64) of [§. 
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